Quantum Phase Transition in Finite-Size 
Lipkin-Meshkov-Glick Model 
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Lipkin model of arbitrary particle-number N is studied in terms of exact differential-operator 
representation of spin-operators from which we obtain the low-lying energy spectrum with the in- 
stanton method of quantum tunneling. Our new observation is that the well known quantum phase 
transition can also occur in the finite- N model only if TV is an odd-number. We furthermore demon- 
strate a new type of quantum phase transition characterized by level-crossing which is induced by 
the geometric phase interference and is marvelously periodic with respect to the coupling parame- 
ter. Finally the conventional quantum phase transition is understood intuitively from the tunneling 
formulation in the thermodynamic limit. 
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The Lipkin-Meshkov-Glick (LMG) model of many- 
body two-level system 0, 0, which was invented orig- 
inally for the study of nuclear giant monopole 4] , is ex- 
actly solvable|a, la, [?J and has wide applications in the 
statistical mechanics of mutually interacting spins 0, E3 
and the Bose- Einstein condensates Recently, it has 
been regarded as an important model to explore the 
novel relation between the quantum entanglement and 
quantum phase transition (QPT)[H 111 fll 111 Il5|. 
QPT describing structural change of the ground-state 
energy-spectrum of many-body systems associated with 
the variation of coupling parameters has attracted con- 
siderable attentions in the modern theoretical and ex- 
perimental communities 16] . As a matter of fact QPT 
originates from singularity of the energy spectrum |l7j 
and any non-analyticity point in the ground-state en- 
ergy of the infinite-lattice-system can be identified as a 
QPT. There can be a level-crossing where the excited- 
state and ground-state levels interchange at the critical- 
point, creating non-analyticity of the ground-state energy 
which is a function of the coupling parameter [16j. It has 
been demonstrated that the LMG model described by 
the Hamiltonian |Ts| 
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undergoes a QPT from a deformed to a normal phases at 
the critical-value of coupling parameter h c = A' = NX, 
where the parameter A' is used to avoid the formal 
infinity |la |20| in the thermodynamic limit TV — > oo and 
a a is the Pauli matrices. In the deformed phase when 
h < h c the ground state is degenerate while is non- 
degenerate in the normal phase when h > h c . For our 
convenience we confine the parameters A>0,0<7<1 
and consider the two cases denoted by " ± " signs sepa- 
rately. In spin language the deformed-phase corresponds 
to a long-range magnetic order and may be regarded 
as a ferromagnet, while the transverse-magnetic- field h 
lowering the potential-barrier induces quantum tunnel- 
ing between two generate-states and eventually destroy 



the magnetic-order at the critical point h = h c . A simi- 
lar QPT was indeed observed with the ferromagnetic mo- 
ment vanishing continuously at a quantum critical point 
of the transverse field in the low-lying magnetic excita- 
tions of the insulator LiHoFe4 |l9j . 

The study of scaling properties at the critical point 
for a finite- TV model is of primer interest and has be- 
come a hard topic recently [IH l2fj| . Based on mean- 
field approximation it is found that particle-number- 
dependence of the energy gap between the ground and 
first-excited states behaviors as exp[— TV] in the deformed 
phase|0, |2(j. However QPT from deformed to normal 
phases, which we called the conventional QPT (CQPT), 
actually does not exist for a finite- TV J8j since the quan- 
tum tunneling can remove the degeneracy leading to the 
tunnel splitting ATS even at h = and thus the ground- 
state energy-structure, strictly speaking, does not have 
macroscopic change. The energy gap AE vanishes in the 
thermodynamic limit TV — > oo and the ground state be- 
comes degenerate. The quantum tunneling effect of the 
finite- TV model, which gives rise to a rich structure of low- 
lying energy spectrum has not yet been studied explicitly. 
Based on the tunneling effect obtained in this Letter we 
demonstrate that the CQPT can occur in the finite-TV 
model if TV is an odd-number due to the quenching of 
quantum tunneling induced by quantum-phase interfer- 
ence and the magnetic-order survives from the tunneling 
till the potential-barrier separating the two degenerate- 
states is suppressed sufficiently low by the transverse 
field. 

Moreover we find a new type of QPT characterized 
by the level-crossing which is induced also by quantum- 
phase-interference. The ground and first-excited states 
interchange alternatively with the increase of h resulting 
in periodic QPT. The finite- TV Lipkin model is studied in 
terms of the exact differential-operator representation of 
the spin-operators which gives rise to a potential-barrier 
tunneling Hamiltonian and the low-lying energy spec- 
trum is obtained by means of instanton method. From 
the viewpoint of tunneling the CQPT can be understood 
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intuitively since the potential-barrier height decreases 
with the increase of h and finally disappears at the critical 
field-value. The level-crossing, however, does not appear 
at all in the CQPT. 

We adopt the giant-spin operators defined by Si = 
~Y^h^\/2 (I = x ,y, z ) to convert the TV-particle model 
into an anisotropic giant-spin Hamiltonian 

H = -\{S 2 x ±"fS 2 y )-hS z . (2) 

with the total spin-quantum-number s = N/2, N/2 — 1, 
N/2 — 2 ■ ••, 1/2, or depending on odd- or even-iV and a 
factor 2 absorbed into the parameters A and h. Because 
of the ferromagnetic interaction between spins, A > 0, 
the ground state and the first-excited states always lie 
in the subspace of maximum spin s = N/2. The 
Hamiltonian of " + " sign describes a giant-spin with an 
easy axis-x (lowest energy) and a hard axis-z (highest 
energy). The external magnetic field h is applied along 
the hard axis, which is critical in generating the QPT 
of level-crossing. For the " — " sign case y-axis becomes 
the hard-axis and the magnetic field is then along the 
medium axis. 

We start from stationary Schrodinger equation 
H<&(tp) = EQ(tp) with Hamiltonian of " — " sign 
and a coordinate-frame rotation such that the hard- 
axis and easy-axis are along z and x respectively for 
the sake of convenience, where tp can be visualized 
as the azimuthal angle of the giant-spin vector and 
the generating function plj $(<£>) can be constructed 
in terms of the conventional eigenstates of S z such 
as $(<£>) = J2m=-s C m exp(zmy))/ \/ (s - m)!(s + mjl, 
which obviously satisfies the boundary condition, + 
2tt) = exp(2iris)^(tp). Thus, we have periodic wave 
functions for even-iV and anti-periodic wave functions 
for odd-iV. The explicit form of the spin oper- 
ator acting on the wave function <fr(tp) is seen to 
bepH S x = scostp — sintp-j^, S y — ssin^ + 

cosv?^, and S z — —ijh- Substituting these differen- 
tial operators into the Schrodinger equation we obtain 
{-a(l - 6sm 2 <p)£s - [X(s - ±)sin2^ - hsmtp}± + 
[A(s 2 cos 2 tp + s sin 2 tp) — hs cos tp]}$(tp) — E<fr(tp), where 
a = A(l +7), S = 1/(1 + 7). I n the new vari- 
able x defined as x(<p) = f£ dp'/ \fl — S sin 2 tp' — 
F((p,S), which is the incomplete elliptic integral of 
the first kind with modulus y/S, the trigonometric 
functions sin tp and cos tp become the Jacobian ellip- 
tic functions sn(x) and cn(x) with the same modulus, 
respectively!^. Using the transformation such that 
$[(j>(x]\ = rfn s (a;)exp{-/itan- 1 [cn(a;)/^7]/2A % /7}V'(a;), 
where dn(x) = yl - 5sn 2 (x) is also a Jacobian elliptic 
function, the stationary Schrodinger equation becomes 

[-a^j + V(x)]ip(x) = Eip(x), (3) 
with the scalar potential given by V(x) = [(cn 2 (x) — 



h(s + \)cn{x) + h 2 /ia]/dn 2 (x), ( = Xs(s + 1) - h 2 /Aa. 
The boundary condition of wave function becomes ip(x + 
AK) = exp(2ms)ip(x) with K denoting the complete el- 
liptic integral of the first kind. In the absence of external 
field h = the scalar potential V(x) is symmetric with re- 
spect to the coordinate origin (see Fig.(l), red-line) with 
potential minima located at x± = ±K in one-period cor- 
responding to the original azimuthal angles tp± — 0, tt 
which are two equilibrium orientations of the giant-spin. 
Disregarding the negligibly small quantum tunneling the 
two ground states are degenerate. The external field sup- 
presses the potential barrier and thus increases the tun- 
neling rate which destroys the magnetic order gradually. 
We now calculate the tunneling rate between the two 
degenerate-states. To this end we begin with the tunnel- 
ing induced Feynman propagator through the potential 
barrier 

< x-,P\x+,-P >=< x-\e~ 2l3 "\x + >= J V{x}e~ S 

(4) 

where S — lim^oo C e dr is the Euclidean action 
evaluated along the tunneling trajectory of a pseudopar- 
ticle in the barrier region called instanton with the Eu- 
clidean Lagrangian C e — x 2 /4a + U(x). x = dx/dr de- 
notes the imaginary-time derivative with t = it. The 
quantum tunneling removes the degeneracy of ground- 
states \ip c ±) m two potential wells located at x± respec- 
tively which describe the two equilibrium-orientations of 
the giant-spin and thus may be called the Schrodinger cat 
states. In the two-level approximation we have H |1) = 
Ex \1), H |0) = E a |0), with |0), |1) = (|Vc+) ± \4>c-))/V2 
being the ground and first-excited states. The low- lying 
energy spectrum can be evaluated as Ei$ = E ± AE/2 

with E =< ip c ±\H\ip c ± >, and AE = E x - E = -(< 
ip c +\H\ip c - > + < ip c -\H\ip c+ >). Inserting the com- 
plete set \ip c ±) in the transition amplitude eq.(4) the tun- 
nel splitting can be derived from the path-integral 

AE ~ 7 r / T>{x}e~ S , (5) 

which can be evaluated in terms of the stationary-phase 
approximation that J D{x}e~ s ~ Ie~ Sc where S c is the 
action along the classical trajectory of instanton x c (t) 
i.e. the solution of classical equation of motion SS = 0. 

I = J J){r]}e~'^^^' : denotes the contribution of quan- 
tum fluctuation around the classical trajectory such that 
x(t) = x c (t) + t](t) with r\ being the small fluctuation. 
Working out the path-integral with both the clockwise 
and anti-clockwise under barrier rotations (see Fig(l)) 
and wave functions 4'c±{x±) noticing the boundary con- 
dition we obtain the level splitting 

AE = Pe-WyJ cosh \ + cos 2ns (6) 

with P = 2 7 / 2 {U%\{l + 7 )(1 - g 2 ) 5 /[7/(l + 7) + 
g 2 /{\ + 7)] 2 ^ 2 } 1/4 , W = ^o(l+7)AM(l + 
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V(l - g 2 )/(l + 7 )) /(l - V(l g 2 )/( l + 7))], and X = 
2g v /[/ /A7tan - 1 [y/7(l- g 2 )/(l + 7 )/g], where U = 
[Ci - C.3/7 + V(Ci " C3/7) 2 + Cl/71/2 and e = ( 2 /2U 
with Ci - As(s + 1) - ft 2 /4A(l + 7), (2 - h(s + 1/2) 
and C3 = h 2 /4X(l + 7). When li = Owe have \ = 
and the square-root factor in eq.(6), which is a direct re- 
sult of quantum phase interference of wave functions in 
eq.(5) between two tunnel paths (Fig.(l)), is v2 for even- 
TV. The quantum tunneling destroys the magnetic order 
and the ground state |0) becomes paramagnet. There 
is no QPT for finite-TV. However the situation is dif- 
ferent for the odd- TV case where the square-root factor 
in eq.(6) becomes zero leading to vanishing tunnel split- 
ting. The degeneracy of ground states at ft = can- 
not be removed by tunneling similar to the Kramer's de- 
generacy in spin system and the ferromagnetic ground 
states \ip c ±) remain. Increase of the magnetic field sup- 
presses the potential barrier seen from Fig.(l) and en- 
hances the tunneling rate which destroys the magnetic 
order gradually. The system thus undergoes a smooth 
second order QPT. The critical transition point h c = 
A(l - 7 )(s + 1/2) + \y/(l - 7 ) 2 (s + 1/2) 2 + 4 7 s(s + 1) 
can be determined precisely from the condition that the 
central barrier-height vanishes and the double-well of the 
potential becomes a single one (see Fig(l)), which reduces 
exactly to ft*' 1 = A' in the thermodynamic limit [Tsl. f20j| . 
The numerical simulation of the low-lying levels is shown 
in Fig. (2) the energy gap obtained from which agrees 
with our gap formula eq.(6) perfectly for any s. From 
eq.(6) the finite-size scaling behavior of the energy gap 
is AE ~ exp[— TV] which vanishes in the thermodynamic 
limit no matter TV is even or odd. 

We now consider the most interesting case of " + " sign 
where the external field is along the hard-axis (z-axis). 
The Hamiltonian operator does not possess a satisfac- 
tory differential-operator representation[23|. The Feyn- 
man propagator eq.(4) in angle variables may be evalu- 
ated with the spin-coherent-state path-integrals with the 
spin coherent state defined by 5-n|n >= s|n > where 
n = ( sin 9 cos tp, sin 9 sin tp, cos 9) is a unit vector [^J [25| . 
Regarding tp, p — s cos 9 as canonical variables and 
integrating over variable p we have the effective La- 
grangian given by C = \ra{tp)tp 2 + A(tp)tp — V(<p), where 
m(tp) = 1/2A(1 - (1 - 7) sin 2 tp), A(<p) = s[l - ft/2As(l - 
(I-7) sin 2 tp)],imdV(<p) = \{l-i)s{s+l) sin 2 p-h 2 {\- 
7) sin 2 ip/4X(l — (1 — 7) sin 2 tp). The periodic potential 
V(tp) has degenerate vacua located at tp = 0, n cor- 
responding to the equilibrium orientations of the giant 
spin along ±x directions. We again consider the tunnel- 
ing induced propagator between two degenerate vacua. 
The total-timc-dcrivativc term A(tp)tp coming from the 
geometric phase of the path-integral does not affect the 
classical equation of motion but gives rise to the quan- 
tum phase interference between the clockwise and anti- 
clockwise tunneling paths. With the same procedure the 
tunnel splitting is found as 

AE = Pe- W cos[tt(s - 9 h )] (7) 



with P = 2Vs 3 A 2 (l - 7) 3 / 2 (l - ?7) 5 / 2 /7r(7 - 77), 



W = y/s(s + l){\n{(VT=7j + y/T^)/(^T=7j - 

vT=7)] - y/vfiMO- + vW(l - v^)]} and Oh = 
h/2\yfy, where r\ = h 2 /A\ 2 s(s+l) and ? = (1-7)17/7(1- 
77). We now come to the most interesting observation 
of the paper that the energy gap is a periodic function 
of h and the two states interchange periodically. The 
the QPT of level-crossing arises at the critical points 
h c (n) = 2An % /7 for odd-iV, (2n + l)A^/7 for even-TV, 
where n is integer. Fig. (3) shows the numerical result 
of the first-excited (red-line) and ground (blue-line) en- 
ergy levels for odd-A^ (a) and even-TV (b) respectively. 
The absolute value of numerical energy-gap as a func- 
tion of h agrees qualitatively with eq.(7) and the agree- 
ment becomes perfect for large TV (see Fig(3)-insert) since 
the spin-coherent-state-path-integral method is valid for 
large s. For the odd-TV case the ground state at ft = 
(6h = 0) is degenerate since cos(7rs) = 0. When 
h < ft c (l) quantum tunneling destroys the degeneracy 
and ground state becomes quantum paramagnetic while 
the first excited state is ferromagnetic. This can be ex- 
plained with the Schrodinger cat states that \ip c ±) = | n ±) 
where S , 2; |n± >= ±s|n± > with the unit vector of spin- 
coherent-state along the x-axis. We have < 0|5a;|0 >= 0, 
and < llS^l >= s. Increase of the coupling h cross- 
ing over the first critical point h c (l) the first-excited and 
ground states interchange and system becomes the fer- 
romagnetic phase. The interchange takes place again 
at the second critical point and is particularly useful 
in the quantum computing since the two phases can be 
used to realize a qubit in some practical spin-system. 
The magnetization as a function of the magnetic field 
can be evaluated as Mo,i = — < 0, l\dH/dh\0, 1 >= 
—dE i/dh. The numerical simulation of the magne- 
tization is shown in Fig. 4. The critical behavior of 
the energy gap AE at the critical point h c is found 
as AE(h -► h c ) = {Tr/X^i)P{h c )e~ w{h ^ \h - h c \ . In 
the thermodynamic limit the tunnel splitting vanishes 
as linijv^oo AE ~ exp[-TV] -> and the QPT of level- 
crossing disappears. A critical value ft*' 1 = A' a/7 of CQPT 
can be determined again from the vanishing barrier- 
height, which differs from mean-field value by a factor 

V7- 

In conclusion: with the explicit calculation of tunnel- 
ing effect we show that QPT can occur in finite-size LMG 
model. The QPT of level-crossing discovered in this pa- 
per may provide a possible technique to realize qubit in 
quantum computer and a paper explaining this procedure 
will be published elsewhere. 
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Figure Captions 

Fig. 1 The potential barriers and clockwise and anti- 
clockwise tunnel paths for h — (red line) and 5 with 
7 = 0.5 in the unit A = 1. 

Fig. 2 The numerical simulation of the first-excited (red- 
line) and ground-state (blue line) energy levels as a func- 
tion of h for N = 11 and 10 (insert) with 7 = 0.95. 
Fig. 3 The numerical simulation of the first-excited and 
ground-state energy levels as a function of h for N = 11 
(a) and N = 10 (b). Insert: the analytical (solid line) and 
numerical (dotted line) absolute values of energy gap for 
N = 41 (a) and N = 40 (b). 

Fig. 4 The first-excited- (red line) and ground-state (blue 
line) magnetizations M as a function of h with 7 = 0.95, 
N = 10. 



